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ABSTRACT 


This  report  susmarizes  work  performed  under  Contract  H00014-70-A-0166, 
Task  0001,  Item  9,  Exhibit  I,  Sonar  Array  Studies.  The  study 
described  was  aimed  at  developing  design  guidelines  for  sonar 
projecting  arrays  to  be  used  l'or  high  source  level ,  wide  sector 
applications.  The  effort  began  with  a  literature  survey  described 
in  the  first  two  sections:  "fieneral  Power  Limits  of  Sonar  Trans¬ 
ducers,"  and  "Piezoelectric  Properties  of  Ceramics."  The  piezoelec¬ 
tric  equations  of  state  are  presented,  along  with  a  discussion 
of  and  results  from  two  computer  programs  written  to  predict  the 
performance  of  some  different  ceramics  under  various  environmental 
and  operational  conditions.  An  analysis  of  beam  patterns  and 
directivity  index  of  curved  face  arrays  is  presented  in  Section  V. 

The  final  section  contains  conclusions  and  recommendations. 
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HEFERHMCES 


I.  INTRODUCTION 


This  re?-"". t  tiar.orizes  the  work  performed  under  Contract 
N00014-70-A-0166,  Task  0001,  Item  9,  Exhibit  I,  Sonar  Array  Studies. 

The  purpose  ol'  this  work  was  to  investigate  curved  face  array  configu¬ 
rations  to  provide  guidelines  for  designing  high  acoustic  power,  high 
search  rate  sonar  transmitting  arrays  capable  of  large  angular  sector 
coverage.  High  search  rate  sonars,  operating  at  high  power  levels, 
increase  the  probability  of  ceramic  element  failure  resulting  from 
overheating  and/or  overstressing.  Such  instances  have  occurred  with 
some  prototype  and  experimental  projectors  built  at  Applied  Research 
Laboratories  (AHL).  If,  as  proposed,  ’  future  generation  transmitting 
arrays  are  to  be  operated  at  even  higher  search  rates  and  with  an 
order  of  magnitude  increase  over  present  operating  j>ower  levels,  it 
seems  imperative  to  know  the  effects  of  these  future  requirements  on 
sonar  performance,  so  that  adequate  design  guidelines  may  be  formulated. 
This  report  documents  the  results  of  a  one  year  effort  directed  toward 
prediction  of  projector  performance. 


To  limit  the  problem  to  a  manageable  scope,  the  present  effort 
was  concerned  only  with  piezoelectric  ceramic  arrays.  The  dismissal  of 
magnetostrictive  sources  is  at  least  partially  Justified  in  an  effort 
to  obtain  maximum  acoustic  output  power: 


For  low  Q  transducers,  heavily  loaded  by  radiation  (low  Q  ) 
the  presently  available  nickel-cobalt  alloys  are  capable  of 
far  less  power/pound  than  are  the  lead  zirconate  titanate 
ceramics.  On  the  other  hand,  the  strength  of  the  magneto¬ 
strictive  metals  is  an  important  advantage,  especially  for 
large  structures  such  as  low  frequency  ring  transducers  .... 
Breakthroughs  in  fundamental  research  would  probably  be 
required,  however,  to  put  magnetostrictive  materials  in 
a  position  of  general  supremacy.? 
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Although  hydroacoustic  sources  of  sound  are  suited  for  low 
frequency  applications,  the  hydroacoustic  projector  was  not  considered 
since  the  primary  interest  in  this  study  was  the  ultrasonic  frequency 
range. 

The  study  was  conducted  as  parallel  efforts  in  two  areas:  array 
configuration,  and  properties  and  limitations  of  the  piezoelectric 
ceramic  materials. 

The  remainder  of  this  report  is  divided  into  five  major  sections. 
Section  II  deals  with  general  power  limitations  of  sonar  transducers. 
Section  III  briefly  reviews  the  important  piezoelectric  properties  of 
ceramics  and  discusses  some  important  factors  which  affect  the  various 
properties. 

In  Section  IV,  a  theoretical  analysis  of  transducer  performance 
based  on  the  piezoelectric  equations  of  state  is  presented.  A  mathe¬ 
matical  model  is  developed  for  a  lossless,  thickness -mode  piezoelectric 
plate,  and  the  equivalent  circuit  technique  is  demonstrated.  The 
analysis  is  extended  to  include  both  mechanical  and  electrical  losses 
in  the  ceramic.  Details  of  two  computer  programs  based  on  this  analysis 
are  discussed. 

An  analysis  of  curved  face  sonar  arrays  is  given  in  Section  V. 
Beamwidth  and  directivity  index  calculations  are  considered  os  a 
function  of  array  configuration,  element  size  and  type,  and  amplitude 
shading  techniques.  A  relationship  between  directivity,  source  level, 
and  acoustic  power  is  also  given.  The  final  section  (VI)  contains 
conclusions  end  recoapendatlons . 
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II.  CENERA1  POWER  LIMITS  OF  SONAR  TRANSDUCERS 

The  power  output  of  a  sonar  transducer  nay  he  considered  to  have 
an  electrical,  a  mechanical,  an  acoustical,  and  a  thermal  limit. 

WoolletP  discusses  the  electrical  and  mechanical  limits  for  both  an 
electric  field  transducer  and  a  magnetic  field  transducer.  The  elec¬ 
trical  limit  is  proportional  to  storable  electric  and  elastic  energies, 
mechanical  Q,  mechanoacoustical  efficiency,  and  electromechanical 
coupling  coefficient.  Woollett's  analysis  was  done  for  linear  electro¬ 
acoustic  transducers.  The  limits  on  stored  energy  are  Imposed  by  such 
factors  as  insulation  breakdown,  depolarization  of  materials,  distortion 
resulting  from  dielectric  (and  ferromagnetic)  nonlinearities,  and 
deterioration  of  efficiency  resulting  from  a  rise  in  the  dissipation 
factor  at  high  fields.  Thus,  it  would  seem  that  a  dynamic  approach  to 
power  limits  would  be  required  in  considering  the  high  power  (and  there¬ 
fore  high  driving  field)  operating  conditions.  The  mechanical  limit 
depends  on  the  transducer's  capacity  to  store  elastic  energy  and  is 
primarily  determined  by  the  elasticity  of  ceramic  or  crystal  transducers. 

UrlcP  discusses  the  cavitation  limitation  of  acoustic  power.  The 
expression  relating  cavitation  threshold  and  sound  pressure  level  is 


1(0.707) 

X  =  J - 


106  ?  7  P  ? 

: - Si-  X  10''  =  0.3  Pc  W/cm  , 


I  =  cavitation  threshold,  V/cm  , 

=  peak  pressure  of  sound  wave  causing  cavitation,  atm, 
p  =  density  of  water,  g/cm^, 
c  =*  sound  speed  in  water,  m/s ec. 


The  cavitation  limit  for  a  given  projector  can  be  obtained  by 
multiplying  Ic  by  the  surface  area  of  the  transducer  and  using 
adjusted  for  operating  depth.  The  cavitation  threshold 


Ic  =  0.3  Pe‘  A  , 

where  A  is  the  surface  area  of  the  projector  (cm  ),  represents  the 
maximum  working  value  for  the  power  output  (W)  of  the  projector.  At 
the  limiting  power  level  several  phenomena  begin  to  occur  which 
degrade  the  projector  performance.  These  include  erosion  of  the  pro¬ 
jector  surface,  loss  of  acoustic  power  In  absorption  and  acattering 
by  the  cavitation  bubble  cloud,  a  deterioration  in  the  beam  pattern  of 
the  projector,  and  a  reduction  of  the  acoustic  impedance  into  which  the 
projector  must  operate.  The  cavitation  threshold  is  inversely  propor¬ 
tional  to  pulselength  and  directly  proportional  to  frequency.  An 
approximation  to  measured  values  is  shown  in  Dwg.  AS-71-668. 

The  sound  pressure  at  the  cavitation  threshold  is  a  function  of 
depth  (hydrostatic  pressure), 

Pc(h)  =  Pc(0)  +  ~  atm 

where  h  is  in  feet  (33  ft  of  water  is  equivalent  to  1  atm  of  pressure). 
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Sherman  introduced  a  factor  7  to  account  for  the  nearfield 
effect  on  the  cavitation  threshold, 


i0  =  0.3  7(pc«>)  ♦  y2  . 


Experimental  values  for  y  are  betveen  O.J  and  0.6  for  a  circular  piston, 
vith  an  average  of  about  0.5.  Therefore,  the  cavitation  threshold  of 
a  transducer  can  he  expressed  as 

Ic(f,h,T,A)  =  0.15^0)  +  .  g(f)  .  e(T) 

where  g(f)  is  frequency  dependence  and  e(r)  is  the  pulselength  factor, 
as  described  in  Dwg.  AS-71-668. 

The  thermal  limit  of  acoustic  power  is  exceeded,  for  a  transducer 
composed  of  piezoelectric  ceramic,  when  the  power  dissipated  within 
t.he  element  is  sufficient  to  raise  the  temperature  of  the  material  to 
the  Curie  temperature.  At  the  Curie  temperature,  the  piezoelectric 
property  is  destroyed  and  is  not  restored  when  the  ceramic  cools. 

The  temperature  behavior  of  a  piece  of  piezoelectric  ceramic  is 
a  function  of  time,  of  the  thermal  conductivity  of  the  surrounding  media, 
of  ceramic  volume,  of  ambient  temperature  signatures,  and  of  the 
dissipation  factor  of  the  ceramic.  The  dissipation  factor  (tan  8)  is 
the  ratio  of  energy  dissipated  within  the  ceramic  in  the  form  of  heat 
to  we  energy  converted  into  mechanical  motion.  It  varies  with  planar 
stress  (hydrostatic  pressure),  driving  field,  and  type  of  ceramic. 

Only  a  limited  amount  of  experimental  or  theoretical  work  has 

been  published  on  the  thermal  behavior  of  electroacoustic  projectors, 

g 

Smith  and  Burbage  performed  some  experiments  in  connection  with  a 
comparative  study  of  piezoelectric  ceramics.  They  found  that  specimens 
of  commercially  available  ceramics  very  similar  in  low  field  properties 
can  show  extreme  variation  in  high  field  loss  behavior.  One  anomaly 
observed  at  high  operating  fields  was  the  occurrence  of  "hot  spots," 
uneven  heating  of  the  ceramic. 
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HI.  PIEZOELECTRIC  PROPERTIES  OF  CERAMICS 


The  power  level  capability  of  a  projector  is  dependent  upon 
properties  of  the  piezoelectric  ceramic.  The  values  of  these  parame¬ 
ters  are  available  from  data  sheets  published  by  various  suppliers . 
Typical  values  for  piezoelectric  ceramics  conaonly  used  in  projectors 
ire  given  in  Table  I. 

Several  studies  of  ceramic  behavior  under  various  operational  and 
environmental  conditions  have  been  published  in  the  literature. 

Drawings  AS-71-669  through  AS-71-6f4  are  adaptations  of  experimental 
data.  The  ceramics  are  referred  to  by  brand  name,  and  only  the  behavior 
of  projector  materials  is  reported  here. 

9 

Brown  tested  the  variations  of  several  properties  of  PZT4  and 
NRE4  ceramics  as  a  function  of  planar  stress  and  electric  field. 

Drawings  AS-71-669,  AS-71-670,  and  AS-71-671  are  adaptations  of  his 
results,  using  one-dimensional  static  pressure  for  conventional  mass 
loaded  stacks.  Martin10  reported  on  the  variation  of  coupling  factor 
with  pressure  and  his  results  are  shown  with  Brown's  in  Dug.  AE-71-672. 

Berltneourt11  reported  the  behavior  of  dissipation  factor  as  a 
function  of  electric  field  and  temperature  fer  the  three  Clevlte 
projector  materials,  rZTU,  FZT8,  and  Ceramic  B.  These  results  have 
been  suasarized  in  Dug.  AS-71-671, 

12 

Gerson  observed  variations  of  mechanical  Q  and  Young's  Modulus 
as  a  function  of  stress  (Dug.  AS-71-674).  Hie  work  was  done  with  zvo 
ceramic  compositions  with  characteristics  different  from  the  so-called 
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Y  OF  PROJECTOR  CERAMIC  PROPERTIES 


RELATIVE  PIEZOELECTRIC  RELA f IVE  DIELECTRIC 

CONSTANT  d/dn  CONSTANT  K/K0 


J - 1 - 1 

1  2  3 

PLANAR  STRESS  (x  10*  N/w3) 


VARIATION  OF  DISSIPATION  FACTOR  KITH 
PLANAR  STRESS  AND  ELECTRIC  FIELD 
(NRE4  ond  P2T4)  (9) 


ARL-UT 
AS  .71 -S70 
JMC-M 
J-SS-T1 
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FRACTION  OF  INITIAL 


VARIATION  OF  PERMITTIVITY  WITH  PLANAR  STRESS 
AND  ELECTRIC  FIELD  FOR  NRE4  AND  PZT4  CERAMIC  (9) 


AUL  -  UT 
M-71-67I 
JHC-OR 
5M-7I 
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AAL-UT 
A* -71. *7} 


VARIATION  OF  COUPLING  FACTOR  WITH  PRESSURE  FOR  NRE4 
AND  PZT4  AS  MEASURED  BY  BROWN  (9)  AND  MARTIN  (10) ' 


projector  materials.  No  such  data  have  been  found  on  this  particular 
type  of  ceramic  used  in  projectors.  Cerson's  results  are  reproduced 
here  only  as  an  example  of  what  may  happen  under  certain  conditions. 


XV.  TRANSDUCER  ANALYSIS  EASED  ON  PUZOELECTRIC  EQUATIONS  OF  STATE 


i 
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In  this  section,  the  piezoelectric  equations  of  state  will  be 
used  to  develop  a  mathematical  model  of  a  thickness  code  piezoelectric 
plate  in  which  both  electrical  and  mechanical  losses  are  taken  into 
account.  Once  the  equations  have  been  developed  for  the  general  case, 
an  equivalent  circuit  will  be  developed  for  the  Idealized  case  of  a 
thickness  mode  ceramic  with  no  internal  losses.  Several  examples  are 
cited  from  the  literature  in  which  equivalent  circuit  techniques  were 
used  to  obtain  design  guides  for  transducer  performance  (resonance 
frequency,  power  input,  power  output,  dissipated  power,  impedance).  In 
the  remainder  of  this  section,  two  computer  programs  written  to  predict 
transducer  performance  are  discussed.  The  programs  Involve  an  exact 
solution  to  the  piezoelectric  equation*'  developed  at  the  beginning  of 
this  section.  The  first  program,  CERAMIC,  predicts  transducer  perform¬ 
ance  as  a  function  of  frequency  and  drive  level.  The  second  program, 

C ALTOS',  predicts  the  temperature  rise  and  associated  effects  on 
ceramic  performance  resulting  fraa  the  calculated  power  dissipated  in 
the  ceramic. 


A.  3-Port  Equations  for  Thickness  Mode  Piezoelectric  Plate  with 
Losses  Included  ~ 


A  cannon  mode  for  driving  sonar  projectors  is  that  in  which  the 
applied  electric  field  is  parallel  to  the  direction  of  elastic  wave 
propagation,  and  the  crystal  is  expanded  in  the  thickness  direction 
as  shewn  in  Dwg.  AS -71-1101. 

Pnceiiig  pigs  link 
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It  Is  assumed  (l)  that  the  length  and  width  dimensions  are 
sufficiently  large  so  that  displacement  occur  only  in  the  thickness 
direction,  (?)  that  the  two  plane  surfaces  of  area  1  xv  are  equj- 
potential  surfaces,  and  (3)  that  the  magnitudes  of  applied  or  produced 
electrical  or  mechanical  stresses  are  small  enough  to  insure  linear 
operation. 

I?  mechanical  losses  proportional  to  the  strain  rate  are  assumed, 
then  the  equation  for  the  generated  stress  may  be  written 

T  =  cD(l+a  tan  B  )S-hD  , 


where 


c  is  the  stiffness  constant, 

S  is  the  strain, 
h  is  a  piezoelectric  constant, 

D  is  the  electric  displacement, 

J  is  -/3l,  and 

tan  B  is  the  mechanical  "loss  tangent"  of  the  ceramic, 
c 


By  definition,  the  loss  tangent  is  the  ratio  of  the  elastic  energy 


where  €*  is  the  ccraplex  dielectric  permittivity  defined  as 
55 

=  Sjj  {1-j  tan  6g) 


She  quantity  is  the  real  permittivity  measured  at  constant 
strain  and  tan  8^  is  the  electrical  loss  tangent  for  the  ceramic. 
The  (peak)  applied  voltage  is  given  by 


Use  face  velocities  and  the  input  current  obtained  for  the  3-port 
equations  are  (time  dependence  is  suppressed): 


Z  -  .<  _L_ 

33  3  oe^jA  • 

Detailed  derivations  of  these  equations  appear  in  contract  status 
reports  for  1  December  1970  -  30  May  19711  and  1  June  1971  - 
31  August  1971. 15 

The  3-port  equations  are,  within  the  Units  of  the  assumptions 
nade,  an  exact  description  of  the  kinetic  energy,  the  potential  energy, 
and  the  motion  in  the  mechanical  vibrating  system  as  a  function  of 
frequency.  Because  most  engineers  are  more  familiar  with  electrical 
systems  than  with  mechanical  systems,  mechanical  systems  are  often 
analyzed  in  terms  of  an  equivalent  electrical  circuit.  An  electrical 
circuit  is  analogous  to  a  mechanical  system  if  the  differential  equation 
of  the  circuit  is  the  same  type  as  the  differential  equation  which 
describes  the  mechanical  system  and  if  the  two  systems  exhibit  Identical 
energy  characteristics  at  each  frequency. 

B.  Equivalent  Circuit  of  Ideal  (Lossless)  Piezoelectric  Plate 

To  illustrate  the  application  of  equivalent  circuit  techniques 
and  the  way  in  which  such  analyses  can  yield  useful  design  information, 
the  equivalent  circuit  for  a  thickness  mode  plate  with  no  losses  is 
developed.  Also,  some  results  obtained  through  equivalent  circuit 
analysis  by  various  authors  are  reproduced  in  this  section. 

The  equivalent  electrical  circuit  is  obtained  by  assuming  that 
voltage  is  analogous  to  force,  current  to  velocity,  mass  to  inductance, 
capacitance  to  mechanical  compliance,  and  resistance  to  mechanical 
resistance  (the  so-called  "impedance"  or  "classical"  analogy).  Under 
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these  assumptions,  the  equivalent  circuit  for  the  lossless,  thickness 
node  plate  can  be  shovn  to  be  the  circuit  given  in  Deg.  AS-71-1248, 
where 


C  .?& 
co  l 


and 


0  «  hjjCo  is  the  ideal  transformer  voltage  ratio  (the  lapedance 
ratio  is  02). 


EXACT  EQUIVALENT  CIRCUIT  FOR 
THICKNESS  MODE  PIEZOELECTRIC  PUTE 


ARL, .  UT 
AS  -71  -IMS 
LFH-ORS 
12  *  9  *71 


The  circuit  shovn  in  Dwg.  AS-71-1248  is  an  exact  representation 
of  the  piezoelectric  ceramic  operating  in  the  thickness  mode  and  it  is 
valid  for  any  frequency.  In  using  this  circuit  to  describe  a  general 
situation  in  which  the  ceramic  faces  in.**  mechanical  loads,  it  is 
necessary  to  connect  to  the  mechanical  terminals  of  the  circuit  the 
appropriate  mechanical  load  impedances  presented  by  the  loads.  The 
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final  font  for  the  distributed  circuit  with  loads  attached  is  shown  in 
r>wg.  AS-71-1025. 


I  .  FINAL  FORM  FOR  DISTRIBUTED  EQUIVALENT 

I  CIRCUIT  FOR  THICKNESS  MODE  PIEZOELECTRIC  PLATE  AS-71-1025 

jj  If  the  resonator  is  operated  near  the  mechanical  resonance 

[S  frequency  tQ  «  V J21,  then  the  distributed  impedances  nay  be  approxi- 

mated  by  lumped  impedances.  It  can  be  shavn^  that  the  tern 
jSZ  cot  PI/2  is  approximated  by  a  series -resonant  circuit  with  a 
[I  compliance  C  given  by 

r  S  ® 


The  jtZQ  tan  pi/2  Is  approximated  by  a  shunt-resonant  circuit  for  which 


C 

s 
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x 
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a 


and 


,  2  2 

4,1  C. 
O  8 


The  equivalent  circuit  near  resonance  is  shown  in  Dwg.  AS -71-1026. 


LUMPED  ELEMENT  EQUIVALENT  CIRCUIT  TOR 
THICKNESS  MODE  PIEZOELECTRIC  NEAR  RESONANCE  M 

C.  Optimst  Design  Considerations  Based  on  Equivalent  Circuits 

The  application  of  electrical  circuit  theory  as  an  aid  in 
analyzing  elect ro-aechmo -acoustical  vibrating  systems  has  been  well 
documented  in  the  literature.  The  representation  of  a  piezoelectric 
crystal  as  an  equivalent  electrical  circuit  element  was  originally 
given  by  Van  Dyke.1®  The  analysis  of  the  piezoelectric  transducer 
driving  a  mechanical  and/or  acoustical  load  has  been  discussed  by 
Nason, Olson, ^  and  Katz,20  as  well  as  others.  Work  of  a  more 


26 


21 

specific  nature  was  conducted  by  Roth,  in  which  the  equivalent  circuit 
of  a  thickness  node  piezoelectric  transducer  was  used  to  determine 

the  electrical  driving  point  iapedance  as  a  function  of  loading  and 

22 

frequency.  Thurston  has  perforated  an  analysis  of  the  transducer  loss 
and  bandwidth  as  a  function  of  aechanical  and  electrical  loading  based 
on  the  equivalent  circuit  for  a  piezoelectric  transducer.  Kossoff2^ 
has  made  use  of  equivalent  circuit  techniques  in  analyzing  the  effects 
of  backing  and  matching  on  the  performance  of  PZT7A  transducer. 

Further  references  '  ’  1  of  a  more  general  nature  are  listed  in 

the  bibliography. 

The  effects  of  electrical  and  aechanical  resistances  on  transducer 
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loss  and  bandwidth  have  been  analyzed  by  Thurston  for  a  transducer 
operating  near  resonance.  Three  cases  are  considered:  (l)  a  trans¬ 
ducer  electrically  tuned  by  a  shunt  coil,  (2)  a  transducer  tumd  by 
a  series  coil,  and  (3)  an  untuned  transducer.  Sane  results  of  his 
analysis,  given  below,  are  valid  for  the  case  of  a  transducer  free  on 
one  end  or  a  transducer  driving  a  sysnetrical  load  on  each  end. 

The  transducer  loss  at  resonance  for  a  shunt-tuned  transducer, 
as  determined  by  Thurston,  is  shown  in  Dwg.  AS-71-681  as  a. function  of 
load  and  source  resistance.  Transducer  loss  is  defined  by 


TL  =  10  log  ■ 


where 

P  »  power  actually  absorbed  in  load  R^, 

p 

max  «  power  which  would  be  absorbed  if  load  and  generator  were 
matched. 

The  source  Impedance  it  the  quantity  Rg.  The  transducer  lots  as  a 
function  of  frequency  is  shown  in  Dwg.  AS -71-684. 
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For  the  case  of  the  series-coil-tuned  transducer,  the  results 
are  essentially  the  same  as  the  shunt-coil  case,  as  long  as  the 
fractional  bandwidth  of  the  transducer  is  small.  The  fractional  band¬ 
width  (FSW)  is  given  by 


where  fg  and  f^  are  the  high  and  the  low  cutoff  frequencies 
respectively. 

For  the  case  of  an  untuned  transducer,  the  transducer  loss  and 
the  approximate  FBW  are  given  in  Dwg.  AS-71-683  as  a  function  of  the 
parameter  a  defined  by 

a  “  • 

The  effects  of  backing  and  matching  on  the  performance  of 
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piezoelectric  transducers  have  been  analyzed  in  a  paper  by  Kossoff. 

The  analysis  is  valid  for  a  thickness  mode  transducer  for  which  the 
loading  and  backing  media  are  large  enough  in  extent  so  that  no  energy 
is  reflected  to  the  transducer.  For  the  nonaatched  cases,  the  trans¬ 
ducer  is  assumed  to  drive  directly  into  the  loading  and  backing  media. 

Computer  calculations  were  carried  out  by  Kossoff  for  a  KT7A 
transducer,  2.5  cm  dlam,  with  a  series- resonance  frequency  at  2  MHz. 
The  transmitting  sensitivity  curves  for  several  conditions  of  hacking 
and  matching  are  presented  in  Dwg.  AS -71-682.  The  transducer  is 
assumed  to  be  driving  into  a  water  load.  For  those  curves  in  which 
the  transducer  was  matched  either  to  the  backing,  the  load,  or  both, 
the  matching  layer  was  assumed  to  be  a  quarter-wavelength  thick  at  the 
transmitting  resonance  frequency, 


v  -  fo 4-  -  Terr  • 


JO 


(f  is  the  series  resonance  of  the  t  and  the  C  circuit  elements  of 
o  mm 

Dvg.  AS-71-1026  (see  p.  26).) 


The  end  "k"  quantities  in  the  figures  are  defined  by  Kossoff 
as  follows : 


*b 

*b 

zmb 


mq 


characteristic  mechanical  impedance  of  baching  medium, 

Vv 

characteristic  mechanical  impedance  of  quarter-wave  section 

used  to  match  transducer  to  backing  medium, 

characteristic  mechanical  impedance  of  loading  medium 

z  /z  , 
q'  mq' 

characteristic  mechanical  impedance  of  quarter-wave  load 
matching  section. 


From  these  analyses  it  is  apparent  that  optimum  power  transfer  occurs 
near  the  resonance  frequency  of  the  transducer  with  the  load  and  source 
impedance  matched.  The  bandwidth  of  the  transducer  is  increased  by 
matching  either  to  the  backing  or  to  the  load,  but  only  at  the  expense 
of  reduced  output  power.  However,  matching  both  to  load  and  to  backing 
produces  a  wide  bandwidth  as  well  as  a  satisfactory  transmitting 
response. 


/ 

D.  Transducer  Performance  Predictions  Based  on  Computer  Solutions 
of  3-Port  Equations  — 


Two  FOHTRAH  computer  programs  which  were  written  to  obtain 
theoretical  performance  data  are  discussed  in  this  section. 


1.  Program  CEPAMIC 

Program  CERAMIC  vae  written  to  obtain  a  solution  to  the  }-port 
equations  as  a  function  of  frequency  for  the  piezoelectric  plate.  The 
elastic,  dielectric,  and  piezoelectric  constants  of  the  ceramic(s) 
(Section  III)  provide  the  necessary  input  data.  The  computer  solutions 
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are  obtained  for  ceramics  operated  at  field  levels  for  vhich  input 
data  are  available.  The  individual  ceramic  parameters  must  be  specified 
for  each  applied  field  level.  In  this  study,  three  ceramics  coraaonly 
used  and  commercially  available  as  projector  materials,  FZT6, 

Cbannelite  5400,  and  Channelite  5800,  were  examined. 

The  complex  expressions  for  the  (time-independent)  face  velocities 
and  the  input  current  are  reprinted  for  convenience; 


The  equations  fcr  V^,  V£,  and  I  are  solved  by  CERAMIC  as  functions  of 
frequency  and  field  level  and  she  values  are  used  to  calculate  the  rms 
input  and  output  powers  and  the  rms  dissipated  electrical  and  mechanical 
powers  as  functions  and  frequency  for  each  field  level.  The  rms 
electrical  input  power  fcr  sr.eady  state  current  and  voltage  is  given  by 

P  =  §®12  III  e 

The  input  power  per 


where  8  is  the  phase  angle  between  and  I. 
volume  of  ceramic  per  kilohertz  is  given  by 


*12  ReJV  V 


cm^  •  kHz 
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where 


Re(Y^n)  is  the  total  input  conductance,  and 

♦12  is  expressed  in  volts,  A.i  in  cubic  meters,  and  f  in 
kilohertz.  The  dielectric  dissipated  power  (P^g),  determined  from  the 
rms  loss  current,  is 


‘DE 


10T  sc 


3J  tan  fcE! 


kHz 


where  (O^g/Z )  is  the  electric  field  magnitude  expressed  in  kV/cm, 


The  ceramic  plate  was  arbitrarily  assumed  to  have  a  radiating 
2 

area  of  0.01  m  and  a  thickness  of  0.015  m.  The  values  of  the  density 
and  sound  velocity  for  the  acoustic  medium  (water)  are,  respectively, 
10  kg/m  and  1500  m/sec.  The  remaining  ceramic  parameters  necessary 
for  the  computer  calculations  are  given  in  Table  IX.  The  velocity  of 
sound  in  the  ceramic  was  calculated  from  the  given  values  of  p  and  3^ 


v 

o 


based  on  the  assumption  that  the  stiffness  constant  is  the  inverse  of 
the  elastic  compliance.  The  piezoelectric  constant,  hjy  is  the  product 
of  the  voltage  constant  and  stiffness  constant. 


The  input  power,  output  power,  and  dissipated  powers  were 
calculated  for  each  i eraaic  as  a  function  of  frequency  for  various 
drive  levels.  For  Kl8,  tbc  applied  electric  field  assumed  integral 
values  of  1,  2,  ...,  7  kV/cm.  The  applied  field  values  ranged  from 
1  to  k  kV/cm  for  Channelite  5*iOC,  and  from  1  to  6  kV/cm  for  Channel- 
ite  5800. 
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As  indicated  in  Table  II,  the  values  of  the  electrical  loss 
tangent,  tan  8g,  and  the  values  of  the  electric  permittivity  were 
specified  as  a  function  of  the  applied  field.  The  table  lists  the 
percent  Increase  (from  the  no-field  condition)  in  the  relative  dielec¬ 
tric  constant,  as  a  function  of  applied  field.  The  permittivity 
at  some  field  strength  E  at  constant  strain  is  calculated  from 
Kjj(O),  and  the  coupling  factor,  k^y  by 

e33(B)  =  4  eo(lS+1)(1-kf3)  • 

The  reasonance  frequency  of  a  broad  plate  resonator  is  defined 
for  the  purposes  of  this  report  as  the  frequency  of  maximum  input 
current.  For  thickness  mode  vibrating  plate  projectors  that  have  a 
constant  D  f'  boundary  condition,  it  is  pertinent  to  note  that  the 
reso_ance  (uency  defined  above  does  not  occur  at  the  frequency  for 
which  the  plate  thickness  is  equal  to  >./ 2,  where  \  is  the  wavelength  of 
the  elastic  disturbance  in  the  ceramic.  The  deviation  of  the  actual 

£3 

resonance  frequency  from  the  half-wave  frequency  is  given  by  Kossoff 
on  the  basis  of  equivalent  circuit  analysis. 


where  CH  is  the  motional  compliance  of  the  near-resonance  equivalent 
circuit  ant  l  is  the  negative  compliance 

ce  *  • 

where  <p»hjjCQ  is  the  electromechanical  tronsfcrmer -turns  ratio. 

The  capacitance  Co  is  proportional  to  the  permittivity,  and 
the  actual  resonance  frequency  will  be  reduced  as  e®,  increases  with 
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Increasing  field  strength.  This  is  illustrated  in  Dwg.  A8-71-990, 
in  which  the  ratio  of  the  actual  resonance  frequency  to  the  half-ware 
frequency,  i»  plotted  as  a  function  of  applied  field.  These 

results  were  obtained  from  computer  calculations  for  PZT8,  Channelite  54C3, 
and  Channelite  5600,  As  can  be  seen  from  the  graph,  the  reduction  in 
resonance  frequency  is  greatest  for  Channelite  5U00  at  all  field  levels. 

The  calctilated  dielectric  loss  power  is  shown  as  a  function  of 
applied  field  in  Dwg.  AS-71-991.  The  electric  loss  tangent  was 
assumed  constant  with  frequency  and  the  results  shown  are  independent 
of  any  frequency  considerations.  The  mechanical  power  losses  are 
plotted  as  a  function  of  applied  field  in  Dwg.  AS-71-992.  These  plots 
were  made  from  data  calculated  at  the  respective  resonance  frequency 
for  each  ceramic  at  each  value  of  applied  field  strength  and,  therefore, 
represent  the  maximum  loss  as  a  function  of  applied  field.  The  last 
set  of  curves  shown  in  Dwg.  AS -71-993  gives  the  calculated  electro- 
acoustical  efficiency  q  as  a  function  of  applied  field  for  each  of 
the  three  ceramics.  The  plotted  values  correspond  to  the  values  of 
T)e&  at  the  calculated  resonance  frequency  of  each  respective  ceramic 
at  each  value  of  applied  field. 

2.  Program  CAUHiP 

The  results  frees  program  CERAMIC  were  obtained  00  the 
assumption  that  all  of  the  energy  dissipated  in  the  ceramic  was  trans¬ 
mitted  into  the  surroundings  so  that  no  change  occurred  in  the  ceramic 
temperature.  To  obtain  a  core  realistic  picture  of  transducer 
performance,  it  .as  essential  to  include  a  phenomenological  description 
of  a  temperature  rise  in  the  ceramic  resulting  from  the  dissipated 
internal  energy.  This  seemed  especially  necessary  because  future  so pay# 
will  probably  operate  at  higher  search  rates  and  at  greater  power 
levels  than  do  present  systems. 
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REDUCTION  OF  RESONANCE  FREQUENCY  AS  A  FUNCTION  OF  APPLIED  FIELD 


Applied  Field  —  k  V/  cm  (Peek) 

MECHANICAL  power  loss  at  resonance  as  a  function  of  applied  field 


ELECTROACOUSTICAL  EFFICIENCY  v*  FIELD  STRENGTH  FOR  PZT8,  CHANNELITE 


Time  was  insufficient  for  detailed  development  of  a  temperature 
model;  hence,  what  is  presented  here  may  appear  somewhat  oversimplified; 
however,  this  analysis  is  a  beginning  from  which  a  more  sophisticated 
description  can  evolve. 

If  5  represents  the  total  power  dissipated  in  the  ceramic 
(in  units  of  W/ar* /kHz),  then  at  a  frequency  f  (kHz),  the  power/ cm^ 
dissipated  in  the  ceramic  is 

PD  -  QfJW/cm3) 

The  dissipated  power  per  kilogram  of  mass  is  then 
1000  Qf  Jp  (W/kg) 

where  p  is  the  density  in  g/cm3.  The  specific  heat,  Cp,  (j/kg/deg) 
can  then  be  used  to  obtain  the  rate  of  temperature  rise  in  the  ceramic 


To  simplify  the  analysis,  it  is  assumed  that  no  heat  transfer  occurs 
from  the  ceramic  to  the  surroundings.  After  a  time  t  the  heat  rise  is 

or  =  1000  —2  T 


and  1 

i 

i 

J  =  T  + 
f  o 


where  Tf  and  Tq  are  the  final  and  initial  temperatures,  respectively. 
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The  temperature  model  is  completed  by  making  the  following 
assumptions: 

1)  the  ceramic  operates  under  steady  state  conditions  during 
the  time  interval  t,  and 

2)  no  temperature  change  occurs  during  the  period  between  the 
end  of  one  time  interval  and  the  beginning  of  the  next. 

Because  of  these  restrictions,  the  model  can  be  interpreted  as  a 
cv  system  with  a  parameter  update  time  r  or  as  a  pulse  mode  system 
with  a  pulselength  t. 

In  program  CALTEMP  this  temperature  model  was  used  to 
predict  transducer  performance,  taking  into  account  as  associated 
ceramic  temperature  rise.  The  ceramic  dimensions  were  assumed  to  be 
the  same  as  in  CERAMIC : 

P 

Area  =  0.01  m  , 

Thickness  «  0.015  m 

The  load  impedance  on  ceramic  face  1  was  assumed  to  be  zero,  and 
calculations  were  obtained  for  the  case  of  an  air  load  on  face  2  as 
well  as  for  the  case  of  a  water  load. 

The  only  ceramic  for  which  parameter  temperature  data  were 
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available  was  PZTk.  The  published  data  consist  of  the  measured  value 

T 

of  the  dielectric  constant  at  constant  stress,  K^,  and  the  electric 
dissipation  factor,  tan  5g.  Both  parameters  were  measured  at  several 
temperatures  for  field  levels  of  1,  2,  and  3  kV/cra.  Using  the  published 
data,  each  parameter,  at  each  field  level,  was  estimated  for  each  of 
36  temperatures  (the  set  [Tg]),  and  the  results  were  read  into  the 
program  as  a  3  x  36  array  for  each  parameter.  A  1  x  36  array  of  the 
basic  corresponding  temperature  set  [Tg]  used  in  the  calculation  wos 
also  read  in  and  will  be  referred  to  as  the  array  T  . 


The  relative  dielectric  constant  at  constant  stress  vas 
converted  to  the  value  at  constant  strain  from  the  relationship 


where  is  the  coupling  factor.  The  initial  ceramic  temperature,  Tc, 
was  taken  to  he  25*C  for  each  field  level.  The  resonance  frequency  at 
that  temperature  was  determined  for  each  field  level  from  program 
CERAMIC.  These  frequencies  were  used  in  CAI.TCMF  as  the  operating 
frequency  fQ(t.)  for  the  appropriate  field  level  L. 

In  CALT04P,  the  expressions  for  V^,  Vg,  and  1  are  evaluated 
at  the  frequency  fQ(l>)  for  each  pulse  and  the  input,  output,  and  dissi¬ 
pated  powers  are  calculated.  The  temperature  Tf  at  the  end  of  each 
pulse  is  determined  and  compared  with  the  temperature  values  in  the 
array  T^.  The  parameter  values  used  for  the  next  pulse  are  those 
corresponding  to  the  temperature  in  the  array  Tg  nearest  the  calculated 
ceramic  temperature  T^. 

The  results  obtained  from  CALTUfP  for  each  sequence  of 
500  pulses  (updates)  are  presented  in  Dwgs.  AS-71-1315  through 
AS-71-1320.  The  first  three  sets  of  curves  were  obtained  for  bilateral 
radiation  into  a  zero  impedance  load  with  pulselength*-  (parameter 
iqjdates)  of  100,  300,  and  700  msec.  The  temperature  versus  pulse 
number,  input  resistance  versus  temperature,  input  conductance  versus 
temperature,  and  input  power  versus  temperature  are  plotted  for  each 
value  of  pulselength.  A  second  set  of  similar  graphs  is  shown  for  the 
case  of  air  backing,  with  the  ceramic  radiating  into  a  water  load.  In 
each  set  an  additional  curve  is  shown  depicting  the  output  power  as  a 
function  of  temperature. 

The  temperature  effects  become  more  dramatic  the  longer 
the  drive  time  and  the  greater  the  field  strength.  As  indicated  free) 
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PERFORMANCE  MEASURES  FOR  WATER  LOAD,  TAU  =  300 
FOR  FiELD  LEVELS  OF  1,  2  AND  4-kV/cm 
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PERFORMANCE  MEASURES  FOR  WATER  LOAD,  TAU-700 
FOR  FIELD  LEVELS  OF  1,  2  AND  4-lcV/cm 


ARL  -  UT 
AS71-1317 
LHP  -  8JW 
12.  10-71 
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PERFORMANCE  MEASURES  FOR  AIR  LOAD,  TAU  =  300 
FOR  FIELD  LEVELS  CF  1,  2  AND  4-kV/cn 
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the  resistance  versus  temperature  curves,  there  is  a  reduction  in  the 
input  resistance  as  the  temperature  increases.  The  increase  in 
reactance  with  temperature  indicates  a  shift  in  resonance  to  a  lover 
frequency  as  the  temperature  increases.  Hie  severity  of  these  effects 
on  the  input  power  level  is  seen  from  the  reduction  of  input  power 
as  the  ceramic  temperature  rises. 

Experimental  results  of  a  similar  nature  have  been  obtained 
26 

by  Kendig  and  Clarke.  They  report  a  reduction  of  resonance  frequency 
with  increasing  temperature  and  an  Increase  in  input  resistance  with 
increasing  field  strength.  These  results  are  consistent  with  those 
presented  here. 

However,  Kendig  end  Clarke  also  measured  an  increase  in 
input  resistance  with  increasing  temperature  whereas  CALTEM?  calcula¬ 
tions  would  indicate  just  the  opposite.  In  their  measurements,  Kendig 
and  Clarke  adjusted  their  operating  frequency  with  temperature  so  that 
their  ceramic  element  waa  driven  at  resonance  for  all  measurements. 

In  CAIiTEMP,  the  driving  frequency  was  maintained  constant  at  all  tempera 
tures.  Ho  valid  comparison  of  the  results  for  resistance  versus 
temperature  can,  therefore,  be  made. 
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V.  CURVED  FACE  SONAR  ARRAY  ANALYSIS 

This  section  will  cover  the  teen  patterns  of  various  curved  face 
arrays,  directivity  index  of  arrays,  and  source  level,  as  related  to 
the  array  configuration. 

A.  Wide  Sectors 

Any  beamwidth  can  be  achieved  with  a  planar  array  if  power  levels 
inversely  proportional  to  array  length  {and,  therefore,  amount  of 
ceramic  in  the  projector)  are  acceptable.  One  technique  for  increasing 
the  anount  of  cerasdc  in  a  transducer  without  reducing  the  beamwidth  is 
to  use  a  curved  face  array  configuration.  This  can  be  done  by  either 
using  a  section  of  a  piezoelectric  ceramic  cylinder  or  by  arranging 
thickness  node  staves  vertically  as  though  on  the  surface  of  a  cylinder. 
The  former  is  the  limiting  case  of  the  latter.  As  the  staves  beecoe 
mere  densely  packed,  the  farfield  beam  pattern  calculations  become  the 
same. 


The  farfield  pressure  beam  pattern  of  any  array  of  discrete 
radiators  can  be  expressed  as 


He)  -  i 


E  E(e,k)  •  A*  •  e: 
k»l  * 


,i*(e,k) 
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where 


E (8,k)  ie  the  response  of  the  kth  element  in  the  8  direction, 
is  the  relative  amplitude  of  the  kth  element  and  is 

proportional  to  driving  voltage,  mode  coupling  coeffi¬ 
cient,  and  ceramic  area  of  the  kth  element, 
t(8, k)  is  the  phase  of  the  contribution  from  the  kth  element  in 

the  8  direction.  This  may  he  due  to  electronic  phasing 
end/or  geometry  of  the  array. 

For  the  curved  face  array  of  identical,  discrete,  equally  spaced 
line  sources  driven  in  phase,  the  contribution  of  the  kth  element  in 
the  8  direction  is 


E(6,k) 


Mr  8i°yk) 
r  8in7k 


where 

s  is  the  length  of  each  element,  in  inches, 

X  is  the  wavelength,  in  inches, 
and 

7k  •  e  -  -  *ja  k  »  1,  2,  ...,  n, 

-  8  -  I  +  (k-l)a  , 

with 

N  *  the  number  of  radiators  in  the  array, 
a  =  angular  spacing  between  element  centers 
^2  •  sin*1  aa  shown  in  Dwg.  AS-72-48jj , 
t  o  angular  width  of  the  array  (see  Dwg.  AS-72-l4!j). 

Each  element  la  driven  in  exactly  the  same  manner  so  that 
al,  k®l,  2,  ...,H  • 


5^ 


AS-72->t83 

For  an  array  of  identical  elements  driven  electrically  in  parallel 
under  the  same  coupling  mode  conditions,  the  function  ♦  is  given  by 

*(e,k)  •  ^5-^  cos  rk  , 

where 

\  and  rk  are  aR  defined  above,  and 
R  is  the  radius  of  curvature  of  the  array,  in  inches. 

There  are  several  parameters  vhich  affect  the  behavior  of  F(0). 

In  order  to  isolate  the  effects  of  some  of  them,  the  following  simpli¬ 
fications  were  made: 

(l)  To  eliminate  the  frequency  (wavelength)  dependence,  the 
distance  measures  are  expressed  in  units  of  wavelengths: 

s 

'  “  \  ' 

R’  =  £  ,  and 
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(2)  If  the  element  size  is  allowed  to  approach  zero,  the  — 
term  in  the  animation  is  almost  constant.  In  the  limiting  case, 

lim 
s'-*  0 

for  all  7.  • 

In  the  past,  curved  face  projectors  huilt  at  AHL  have  been 
required  to  insonify  a  maxlimmi  sector  of  about  70*.  As  a  general  rule, 
AHL  curved  face  projectors  have  teen  constructed  with  an  arc  of 
angular  width  approximately  lyf,  to  20 $  larger  than  the  desired  beaawidth 
for  example,  to  achieve  a  60*  beaawidth,  an  arc  of  approximately  70° 
was  used.  This  general  rule  is  also  valid  for  wider  beamwidths.  A 
series  of  patterns  was  calculated  for  different  angular  width  arrays 
and  radii  of  curvature  of  5X,  10\,  SOX,  and  50X.  The  elements  were 
assumed  omnidirectional  to  eliminate  the  effect  of  the  individual 
element  patterns  in  the  calculation.  The  relationship  between  the 
array  parameters  R  and  9  and  the  beaawidth  is  shown  in  Dvg.  AS-71-701. 

With  linear  arrays,  a  zonal  shading  technique  will  broaden  the 
beam  as  well  as  lower  the  side  lobe  levels,  so  that  a  wider  sector  can 
he  covered  without  much  loss  in  directivity  index.  Patterns  utilizing 
several  zonal  schemes  with  curved  face  arrays  were  calculated.  Four 
patterns  are  shown  in  Dug.  AS-71-1086  which  illustrate  the  effects  of 
shading.  All  patterns  are  for  an  array  of  65  point  sources,  equally 
spaced  around  a  130°  arc  with  radius  of  curvature  of  10  wavelengths. 

The  directivity  is  increased  by  these  zonal  shading  techniques,  hut 
this  increase  is  due  to  a  severe  narrowing  of  the  main  beam  rather  than 
side  lobe  suppression. 
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sln(jts*  sin  7^) 
*»'  sin  y. 


x  sin  7j,[  cos(  its 1  sin  7^)11 


IFfECTS  OF  ZONAL  AWirrUOE  SHADING 
ON  THE  BEAN  PATTERN  OF  A  CURVED  FACE  ARRAY 

AS-71-1086 
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In  order  to  determine  the  influence,  if  any,  of  the  type  of 
element  used,  the  same  array  geometry  and  amplitude  shading  were  used 
to  calculate  beam  patterns  using  the  following  individual  radiators: 

(1)  Point 

(a)  Line 

(3)  Piston 

(4)  Circular  Ping. 

It  was  noted  that  for  the  parameters  used  ( cenrer-to-center  spacing 
*  0 . 35\)  the  beam  patterns  do  not  vary  significantly  from  the  patterns 
calculated  with  point  sources.  However,  for  larger  center- to-center 
spacing,  the  point  source  approximation  is  less  accurate . 

In  summary,  then,  any  sector  width  can  be  achieved  using  curved 
face  arrays.  The  limitations  on  the  source  level  are  those  associated 
with  the  ceramic  itself,  available  power,  cavitation  limits,  and  any 
limits  on  the  physical  site  of  the  array.  The  amount  of  ceramic  used 
to  convert  electrical  to  acoustic  power  can  be  increased  by  using  a 
larger  radius  of  curvature  and/or  using  a  longer  vertical  stave  (or 
array),  with  a  correspondingly  narrower  beam. 

A  wider  sector  can  be  covered  in  both  the  horizontal  and  vertical 
planes  hy  curving  the  face  of  the  transducer  in  both  directions. 

Beam  patterns  were  calculated  for  such  a  configuration,  and,  as  one 
might  expect,  the  response  in  the  respective  planes  was  equivalent  to 
those  for  the  cylindrical  configuration  with  the  same  radius  of 
curvature. 

B.  Directivity  Index 

In  order  to  relate  the  beam  pattern  of  a  curved  face  array 
configuration  to  the  source  levels  which  can  be  achieved,  the  relation¬ 
ship  between  the  beamwidth  (or  sector  covered)  and  the  directivity 
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must  be  established.  The  directivity  Index  for  any  transducer  can 
be  calculated  from  Its  three-dimensional  been  pattern.  This  Involves, 
however,  a  numeric  or  analytic  volume  integral  of  the  beam  pattern. 

The  msaerlo  form  is  quite  time  consuming  and  a  general  analytic  expres¬ 
sion  for  the  beam  pattern  of  a  curved  face  array  is  not  readily 
attainable.  Two  approximations  to  the  directivity  index  are  outlined 
below.  Results  are  calculated  over  a  wide  range  of  parameter  values. 

For  wide  angle  curved  face  projectors,  there  are  several  possible 
configurations,  depending  on  the  nature  of  the  desired  pattern  and 
the  amount  of  ceramic  needed  to  handle  the  required  electroacoustical 
power  conversion,  The  two  which  will  be  considered  here  are  the  sector 
of  a  cylinder  (curved  in  one  direction,  planar  in  the  perpendicular 
direction)  and  the  sector  of  a  torus  (i.e.,  two  radii  of  curvature). 

29 

C.  T.  Molloy  7  has  derived  the  directivity  index  for  the  first 
configuration,  using  the  following  simplifying  assumptions: 

1)  The  intensity  of  radiation  is  different  from  zero  only  in 
the  region  -0/2S9S4/2,  ~x/2S(,St/Z. 

2)  Inside  the  above  defined  region,  the  farfield,  the  intensity 
is  assumed  independent  of  9, 

3)  The  dependence  of  the  intensity  on  7  is  assumed  to  be  given 
by  the  following  expression: 

■  ■  hiiuMHA>r  ■ 

where  the  geometry  is  a3smed  to  be  as  shown  in  Dvg.  AS-71-1253. 


The  approximation  can  he  improved  by  using  the  actual  angular 
beunrldth  instead  of  the  sectoral  width  (9)  as  the  boundary  of  the 
region  of  insoaifi cation. 


The  functional  fora  for  directivity  index  (Dl)  is  given  hy 


01  =  -10  lo^  (»(z)  •  a) 

where 

z  *  0.46393  X  (vertical  height  of  staves,  in  inches)  X  f 
(frequency,  in  kilohertz). 

Si  is  the  Integral  sine  function, 

and 

a  is  the  angular  half-power  beamwidth,  in  radians,  in  the  plane 
of  curvature  of  the  horn  or  cylindrical  sector. 

The  values  for  which  the  curves  were  calculated  by  Molloy  were  far  the 
application  of  loudspeakers  in  air.  The  curves  have  been  calculated 
for  values  of  general  interest  in  sonar,  and  are  presented  in 
Dwg.  AS-71-1250. 


The  second  configuration  projector  would  be  more  useful  for  high 
power  applications  in  which  a  relatively  large  angular  sector  is 
desired  in  both  the  horizontal  and  vertical  planes. 

Aasuae  the  intensity  is  uniform  and  maximum  within  a  solid  angle 
defined  by  the  angles  a  and  (5  (as  shown  in  Dwg.AS-71-1251),  and  is 
*ero  elsewhere. 

The  directivity  ratio  (Bp)  is  defined  by  Molloy  as  "...  the 
ratio  of  the  total  acoustic  power  output  of  a  radiator  to  the  acoustic 
power  output  of  a  point  source  producing  the  same  pressure  at  the  same 
point  on  the  axia."^  The  directivity  index  is  given  by 

DI  =  -10  logyjHp  . 
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DIRECTIVITY  FUNCTION  (MOUOY™) 


The  total  acoustic  power  is  the  intensity  integrated  over  uhe  entire 
range  of  angles  a  and  S,  or  the  volume  enclosed  by  the  beam  pattern 
funct.'cn.  The  ratio  of  acoustic  pover  outputs  can  then  be  computed 
from  the  volume  of  the  directional  and  omnidirectional  patterns.  The 


r  ^“^directional  (2oeX3)/a  cp 
P  VoluneJmnid<rectional  (4/3)*  I5  2s 


Dj.  =  -10  log^®,^  =  10  logj^gp  j  • 

This  formula  gives  directivity  index  as  a  function  of  the  product  of 
the  beamwldths  in  the  two  perpendicular  planes  of  the  acoustic  axis 
of  the  transducer.  The  approximation  is  as  good  as  Is  the  assumption 
that  the  volume  of  the  figure  shown  in  Dwg.  AS-71-1251  is  the  volume 
witnin  the  real  beam  pattern.  The  "volume”  of  an  actual  beam  pattern 
is  probably  less  than  the  approximation  bee" use  maximum  source  level 
is  seldom,  if  ever,  achieved  throughout  the  entire  sector. 

Directivity  ii lex  as  a  function  of  the  product  op  is  shown  in 
Dug.  AS-T1-1252.  The  angles  must  be  in  radians  for  use  with  this 
figure. 

C.  Source  l~vel 

The  source  level  of  a  projector  can  be  calculated  from  the 
total  acoustic  power  output  and  she  directivity  index: 

Source  Level  =  71.6  +  10  log  P  +  DI  (dB)  , 
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AS-7T-1252 


DIRECTIVITY  INDEX  AS  A  FUNCTION  OF  SECTOR  COVERAGE 


where 


J 


P  *  total  radiated  acoustic  power,  in  watts, 

DI  »  transmitting  directivity  index  of  the  transducer,  and 
'11.6  is  a  factor  (in  decibels)  to  convert-  radiated  power  in 
watts  to  intensity,  reference  the  intensity  with 
P=1  dyne/  cm 

Drawing  AS-71-1094  shows  the  relationship  between  sound  pressure 
level  and  total  acoustic  power  generated  by  the  transducer,  with 
directivity  index  as  a  paraaeter. 


FUNCTION  OF  ACOUSTIC  POWER  AND  DIRECTIVITY  INDEX 


VI.  CORCLUSIOBS  ATO  HECOMGSTEATIOKS 


A. 


The  acoustic  projection  requirements  -  for  high  power,  high  search 
rate,  vide  sector  sonar  systems  can  be  met  using  curved  face  trans¬ 
mitter  arrays.  If  additional  ceramic  is  needed  to  handle  the  conver¬ 
sion  of  electrical  to  acoustic  power,  a  larger  radius  of  curvature  can 
be  used  to  achieve  comparable  sector  coverage  if  the  arc  length  is 
increased  accordingly.  This  method  of  designing  projectors  is  acceptable 
only  as  long  as  the  physical  size  of  the  transducer  remains  reas enable 
for  the  specific  application. 

A  direct  application  of  planar  array  zonal  amplitude  shading  is 
not  particularly  fruitful  because  the  beanwidth  is  severely  narrowed 
and  because  only  the  highest  anqilitude  elements  can  be  driven  at  full 
power.  If  the  beamwidth  of  a  shaded  array  is  acceptable,  more  of  the 
energy  ia  concentrated  in  the  main  beam  and  the  array  is  therefore 
more  directive. 

In  order  to  achieve  the  sector  coverage  and  source  levels  at 
which  this  study  is  aimed  ( UO*  to  60'  in  either  direction),  signifi¬ 
cantly  more  electrical  power  will  be  required  than  is  used  by  present 
sonars.  If,  for  example,  a  JO'  by  60'  beam  is  to  be  achieved  with  a 
uniform  source  level  of  120  dB,  approximately  10  kW  of  acoustic  output 
power  is  required.  Assuming  50#  energy  conversion  efficiency  (a  typical 
value  for  present  systems),  20  kH  of  electrical  power  would  be  required. 
State  of  the  art  power  smplifiers  can  deliver  such  levels  of  power 


for  duty  cycles  up  to  about  0.5$.*  A  0.5$  duty  cycle  would  provide  a 
system  range  of  800  yd  and  range  resolution  of  25  ft  if  achieved  with 
a  5  msec  pulse  and  pulse  repetition  rate  of  one  per  second.  Getter 
resolution  would  lead  to  lower  duty  cycle  and  even  more  readily 
fulfilled  demands  for  electrical  power.  Higher  repetition  rates  or 
cw  operation  could  greatly  increase  electrical  power  requirements,  but 
the  limitations  of  the  ceramics  examined  here  would  be  encountered 
before  all  available  electrical  power  would  be  used. 

The  adverse  effects  on  ceramic  performance  resulting  from  high 
drive  levels  are  docuaented  in  Section  IV.  In  sunoar-.r,  calculations 
predict  that: 

1)  The  resonance  frequency  of  thickness  mode  ceramics  shifts 
to  lower  frequencies  as  field  level  increased. 

2)  The  electroacouatical  efficiency  drops  with  increasing  field 
levels. 

3)  The  ceramic  tesperature  Increases  with  Increasing  field  level. 

U)  The  increase  in  ceramic  temperature  produces  adverse  effects 

such  as  reduction  of  resonance  frequency  and  loss  of  acoustic  power 
output. 

The  assumptions  made  in  analyzing  the  problm  of  electroacoustic 
energy  conversion  and  associated  source  levels  are  necessary  as  well 
as  instructive  in  a  study  such  as  the  one  described  here.  Its  results 
aid  in  understanding  the  relevant  parameters  and  serve  as  guidelines 
for  future  work. 


Present  work  at  AHL. 


B.  Kecomended  Future  Work 

(1;  Derslop  a  ccKTUter  model  for  transmitting  arrays.  Such  a 
model  would  provide  a  prediction  capability,  not  only  for  the  failure 
of  the  ceramic,  but  also  for  other  aspects  of  array  performance 
quality.  A  coesputer  model  beccaes  especially  attractive  when  its 
cost  is  compared  to  the  cost  of  "trial  and  error"  construction  and 
testing. 

(2)  Conduct  a  theoretical  and  experimental  study  of  heat  transfer 
properties  of  transducer  arrays.  This  study  should  include  the  ceramic 
hacking  material,  the  geanetry  of  the  ceramic,  the  acoustic  window, 
and  the  surrounding  medium.  Whatever  energy  is  not  radiated  as 
acoustic  power  will  be  dissipated  within  the  ceramic  as  generated 
heat  or  given  up  to  mounting  losses.  In  general,  a  rise  in  temperature 
has  two  undesirable  effects:  (c)  it  adversely  affects  ceramic  per¬ 
formance,  and  (t)  if  transmitted  to  the  surrounding  medium,  it  tends 

to  lower  the  cavitation  threshold.  Ho  guidelines  are  known  te  exist 
on  the  thermal  design  of  transducers.  An  initial  effort  in  this  area 
would  in  'ude  a  literature  search  for  thermal  data  on  the  materials  used 
in  various  components  of  transducer  systems . 

(3)  Consider  the  effects  of  cavitation  on  transducer  performance. 

This  Item  Is  important  because  the  onset  o?  cavitation  is  the  ultimate 
limit  of  sonar  source  level.  KLto,  if  the  cavitation  threshold  is 
exceeded,  the  load  into  which  the  ceramic  operates  is  taverely  decreased. 
Such  a  decrees*  causes  an  increase  in  electrical  power  dissipated 
within  the  ceramic,  resulting  in  a  rapid  and  possibly  fatal  rise  in 
te-E'-erature .  !■  simple  model  could  be  developed,  using  a  curve  fitted 

to  published  data.  Such  a  model  would  be  of  limited  scope  and  world 
in  no  way  te  a  complete  formulation  of  the  cavitation  process,  but 
would  be  useful  as  a  genertl  guideline  in  a  study  of  sonar  performance. 
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The  onset  of  cavitation  could  most  readily  he  modeled  as  discrete  and 
complete  at  a  predetermined  pressure  at  the  face  of  the  transducer. 

That  pressure  would  he  dependent  on  frequency,  temperature,  depth 
(ambient  pressure),  and  pulselength.  This  approach  ignores  many 
factors  such  as  particulate  matter  and  dissolved  gases  in  the  water. 

(M  Experimentally  verify  the  model  for  the  transducer  array. 

The  credibility  and  usefulness  of  any  computer  model  is  dependent 
upon  comparison  with  experimental  results.  Two  sources  of  experimental 
data  are  available:  (a)  parameters  of  the  acoustic  projection  phenome¬ 
non  measured  in  the  laboratory  where  strict  control  of  conditions  can 
be  enreised,  and/or  (b)  published  experimental  results. 

(5)  Evaluate  the  usefulness  of  newly  developed  ceramics  with 
a  computer  model.  Two  developments  reported  recently  are  the 

whisker-reinforced  ceramics^1  and  the  ferroelectrics  containing  traces 
32 

of  lanthanum  and  manganese.  The  improvement  in  syrtem  performance 
to  be  realized  by  utilizing  such  materials  could  be  critically 
examined  with  a  computer  model  of  a  projecting  array. 
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